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Abstract 

The random-energy model is studied in the presence of random fields. The problem is solved 

CO . 

exactly both in the microcanonical ensemble, without recourse to the replica method, and in the 
canonical ensemble using the replica formalism. The phase diagrams for bimodal and Gaussian 
random fields are investigated in detail. In contrast to the Gaussian case, the bimodal random 
field may lead to a tricritical point and a first-order transition. An interesting feature of the phase 
diagram is the possibility of a first-order transition from paramagnetic to mixed phase. 
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I. INTRODUCTION 

Spin-glass and random-field models [3( have played prominent roles in the study of 
disordered systems in the last few decades. Although the random-exchange and random- 
field effects are usually considered separately, it has been argued that in proton glasses 
such as Rbi-^NH^H^PC^ 0] it is necessary to take into account the effect of random 
fields generated by the presence of impurities. Another example where the spin-glass and 
random-field effects are present simultaneously is the diluted antiferromagnets Fe a; Zni_ a: F2 

HQ. 

The effect of random fields on the well known Sherrington-Kirkpatrick (SK) model for 
spin glass j?| has been investigated for Gaussian random field Q], bimodal random field [^] 

re low temperature properties of the SK model 

1.2] 

er 

spin-glass model where the effect of random fields can be investigated thoroughly. 
The random-energy model (REM) is probably the simplest spin-glass model 



and trimodal random field |l0|. Since t 
is rather difficult to work out explicitly 



retaining some important properties of the SK model. The REM is related to the gener- 
alization of the SK model to include interaction between every set of p-spins [13]. In the 
p — > oo limit the energies of the spin configurations become independent random variables 
and the model reduces to the REM. 

In this paper we investigate the effect of random fields on the REM. The model is given 
as the p — > oo limit of the Hamiltonian 

T~t — ~ Jii...i p £>h ' ' ' Sip ~ Jo SiSj — 22 HiSi, (1) 

ii<—<i p i<j i 

where Si = ±1 are Ising spins, J^...^ are independent quenched Gaussian random couplings 
with zero mean and variance p\J 2 /2N p ~ 1 , Jo > are ferromagnetic couplings and Hi are 
independent identically distributed quenched random fields. 

The Hamiltonian (0) for p = 2 is the SK model in a random field, whereas for p — > oo 
it reduces to the REM model in a random field. We have solved the problem exactly by 
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two complementary approaches. In section 2 we employ the microcanonical formalism 
to obtain the thermodynamic quantities directly. In section 3 we employ the replica for- 
malism [13] to determine the spin-glass order parameters. In section 4 we study the phase 
diagram for bimodal and Gaussian distribution of random fields. Finally, in section 5 we 
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compare our results with the previous studies on related models and make some concluding 
remarks. 



II. MICROCANONICAL APPROACH 



In this section we solve the model in the microcanonical ensemble Il2l . Let S 
(Si, . . . , Sn) denote one of 2 N spin configurations or the microstates of the system. The 
energy of a given microstate is given by 

E S = H(S) = - Y, ^r-vV" A+ W o(S), (2) 

ii<--<i p 

where H.q denotes the part of the Hamiltonian without random couplings. Since Eg are linear 
combinations of Gaussian random variables Ji v ..i , they are themselves Gaussian random 
variables with mean 

(Eg) = H (S) = E° s , (3) 

and covariance 

J 2 N 



ess' = \(E S - Eg)(E S i - Eg,) 



where 



Qss' 



1 

N 



(4) 



(5) 



is the overlap between the microstates S and 5". In the thermodynamic limit, iV — > oo, the 
energies Es and Eg' of two macroscopically distinguishable microstates S and S' become 
uncorrelated in the p — > oo limit, 

(J 2 N 



p 
Iss' 



0, for p — > oo and \qss'\ < 1- 



(6) 



Thus in the p — > oo limit the energies Eg become independent Gaussian random variables. 
The multivariate probability density is then the product of univariate probability densities 
given by 



f Es (E) 



exp 



(E-E 



0\2 

s) 



NJ 2 



(7) 



V^NJ 2 

Let us consider a given sample, that is, a particular realization of the random couplings 
Ji v ..i . The entropy of the sample is given by 



S(E) = k B lnQ(E), 



(8) 



where 

n(E) =Y,5(E-E S ) 

s 

is the density of states. The average density of states is 

WE)) =J2(5(E - E s )) =Y,f Es (E). 
s s 



(9) 



(10) 



Due to the statistical independence of Es the fluctuations around this average is of order 
(f^-E 1 )) -1 / 2 , and thus completely negligible 12|. 

We can rewrite the average density of states in the form 



(Q(E)) 



dE exp 



NJ 2 



E^o-^). 

s 



We recognize 

n (E)=^S(E -E° s ) 

s 

as the density of states of the system described by the Hamiltonian 7i Q . Therefore 



(n(E)) 



1 



dEo exp 



(E - E ) 2 S (E ) 



NJ 2 



k 



B 



(11) 



(12) 



(13) 



VttNJ 2 J-c 

where 

S Q (E Q ) = k B In n (E ) 

is the entropy of the system characterized by the Hamiltonian Tlo- In the thermodynamic 
limit, iV — > oo, we have 



(14) 



\n(tt(E)) = max 

Eq 



(E — E ) 2 S (Eq) 



Eq is determined by 



1 



NJ 2 k B 
8Sq{Eq) 2k B (E - E 



(15) 



(16) 



Tq(Eq) 8E NJ 2 
where Tq(Eq) is by definition the temperature of the system described by the Hamiltonian 

Ho- 

For energies E such that \n(Q(E)) > the average density of states is very large and the 
fluctuation is negligible. Thus we have with probability 1, 



S(E) = k B In = - kB{E Nj2 E ° )2 + Sq(Eq). 



(17) 



For energies E such that ln(Q(E)) < 0, the average density of states is very small. Thus 
with probability 1 there are no samples with this energy. 



The temperature of the system is given by 

1 dS(E) 2k B (E-E ) 



T(E) dE NJ 2 
which coincides with the temperature of the system described by the Hamiltonian Ho, 

T(E) = T (E ). 

Therefore the energy of the system as a function of temperature is given by 

NJ 2 



(18) 



(19) 



E(T)=E (T) 



2k B T' 



(20) 



where E (T) is the energy of the system characterized by the Hamiltonian 7i Q . The entropy 
as a function of the temperature is 

NJ 2 



S(T) = S (T) 



4k B T 2 ' 

These results are valid above a critical temperature T c determined by 

NJ 2 



S(T C ) = S (T C 



0. 



(21) 



(22) 



4k B T 2 

Below this temperature the system is frozen in its ground state. 

These results are valid for any Hamiltonian Ho. We now particularize for the case where 
the Hamiltonian Ho describes the Ising model with infinite range ferromagnetic interactions 
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Jo 



2N 



(23) 



where in the last passage we have dropped the term Jo/2 that is negligible in the thermo- 
dynamic limit. The quadratic term can be linearized using the identity 



2lY J-oc 

and the partition function will be given by 



^ xe ~\x 2 /2+Xax 



(24) 



s 



lf3J N 



2tt 



/ dm exp < N 

J — oo I 



I I 

■-(5J m 2 + — Y, In 2 cosh /?( J m + F; 



In the thermodynamic limit, N — > oo, the Laplace method gives 



In Zn = N max 



-\(3J m 2 + (In 2 cosh (3( J m + H)) 



(25) 



(26) 



where we have used the law of large numbers to write 



— ^ In 2 cosh/?( J m + #;) = (hi2cosh/?(J m + #)) 



(27) 



where (• • •) denotes the expectation value with respect to the random fields H. Thus the 
free energy is given by 



F = -/TMnZo = N 
where the magnetization m is determined by the equation 



\-J$m 2 - - (In 2 cosh (3(J Q m + H)) 
Z p 



m = (tanh /?( J m + H)) . 



(28) 



(29) 



The internal energy E (T) and the entropy Sq(T) follow from usual thermodynamic relations. 

Applying the general results obtained previously for the system described by the full 
Hamiltonian (0) in the p —>■ oo limit, we obtain for the internal energy 



E 
N 



8 J 2 1 

~ ~ g J o m2 - (H tanh p(J m + H)) 



(30) 



and for the entropy 

s (pjy 



pj m 2 -p(H tanh (3(J m + H) ) + (In 2 cosh (3(J m + H)} . (31) 



Nk B 4 

These results are valid for (3 < (3 C where j3 c is determined by 



5(A 
and 



- p c ((H + J m) tanh (3 C (H + J m)) + (In 2 cosh (3 C (H + J m)) = (32) 



m = (tanh /3 C ( J m + H)) . 
For (3 > (3 C the system is frozen in its ground state. Therefore 



(33) 



E{0) = E(f] c ), S(/3) = 0. 



(34) 



III. REPLICA APPROACH 

In this section we solve the model in the canonical ensemble ^J. We use the replica 
identity for the free energy 

(Z n ) — 1 

-/3F= (In Z) = lim ^ — '- , (35) 



6 



to perform the average over the random couplings J^^...^. To evaluate (Z n ) we introduce 
n replicas of the system a — 1, 2, . . . , n, 

(Z n ) = Tr (e-^ELxWC^)) = Tre"^«, (36) 

where 7i e s denotes the effective Hamiltonian that results after taking the average over ran- 
dom couplings, 



off 



N(pjy 



n 



(37) 



We have dropped terms that vanish in the thermodynamic limit, N — > oo. The nonlinear 
terms can be linearized with the help of the asymptotic relation 



,N\f(a) 



l N\f"{a) 
2tt 



( j iXe NX[f(x)-f(x)(x-a)] 



(38) 



which can be proved for A/" (a) > and N —>■ oo applying the Laplace method. In particular 
for f(x) = x 2 /2 the asymptotic relation reduces to the identity (|24|). Omitting the factors 
that do not contribute to the free energy in the thermodynamic limit, N — > oo, we arrive at 



(Z n ) ~ J IJ dq af3 J l[dm a e-^^\ 



(39) 



where 



a</3 



f = ~\pJ 2 n + \{U\V - 1) E & + U E 



2' v ' ^ R ±ap 2 



a</3 



(40) 



In the N —>■ 00 limit we use the law of large numbers to write the last term as an expectation 
value over the random-field distribution and use Laplace method to evaluate the integral. 
The free energy is then given by the stationary value of the functional 



- = --P-J 2 + lim -{ -PJ 2 (p - 1) Y ql B + -Y\ 



-/T 1 (lnTr exp 



</3 

liPJfp E Cp 1 S a S^ + py(H + J m a )S' 



a</3 



(41) 



where (• • •) denotes the expectation value with respect to the random field H. 
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To compute the free energy we assume 



m n = to, 



(42) 



to be independent of replica indices, and parameterize q a p following the Parisi's i^-step 
replica-symmetry-breaking Ansatz In the n — > limit the free energy functional be- 
comes a function of the magnetization m and the parameters 



<q <qi < ■■■ < qx-i < qi< < 1, 



(43) 



= mo < mi < • - • mx < itlk+i = 1, 



and is given by 



F 

N 



(3J 2 



A 



1 + (p ~ 1) J2( m i+i - m i)?f - PVk 



p-1 



i=0 



, ^0 2 

H m 

2 



- / ^ ( G <7 2 (2/ - -H" - Jom)) #o(y), 

J —CO 



where go(y) is given recursively by 

1 



9i-i(y) 



(3mi 



In 



dy'G a 2(y' - y) exp \Pm9iW)] 



for z = 1,...,K with the initial condition 

0jc(2/) = ^ln(2cosh/5?/). 

G a 2(y) denotes the Gaussian distribution function 

l ( y 2 



GAv) 



where the variances of are given by 



JaV2ir 



exp 



2 J 2 a 2 



(J, 



2 V p-1 



2 Pi p-1 p— 1> 



for 



i = 1, 



(44) 



(45) 



(46) 



(47) 



(48) 



(49) 



We first assume that all the g's are less than one, < go < • • • < Qk-i < Qk < 1- Then 



o 2 — > when p — > oo for z = 0, . . . , iT. Using the expansion 

/■°° / T 2 a 2 d 2 \ J 2 rr 2 

j jy'GAv' - v)f(y') = ™p ( J /(v) = i + ^/"(v) + °(A 



(50) 
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we obtain 



F 
N 



(3J 2 



4 



A' 



K 

E 

i=0 



1 + (P - 1) E( m *+i - m *)^ - Wa 1 



+ ^-m 2 - i(ln2 cosh p(H + J m)) 



-E 

8=0 



1 - (1 — 77ii)(taoh (3(H + J m)) + O(a , . . . , cr^, (T Q a u . . . , a a K 



(51) 



Stationarity of the free energy with respect to the variational parameters gives, in the limit 
p -> oo, 

m = (tanh (3(H + J m)), (52) 

and 

9o = 9i = • • • = Qk = (tanh 2 + J m)). (53) 

Thus we arrived at the replica-symmetric solution where all the q's are identical. The free 
energy in the p — > oo limit is given by 



F 

N 



1 



+ m - ^(ln2cosh/3(# + J m)). 
4 2 



(54) 



The entropy is 

5 {(3Jf 



(3((J m + H) tanh /3 '(J m + H)) + (In 2 cosh (3(J m + if)) . (55) 



Nk B 4 

This solution corresponds precisely to the high temperature solution found in the micro- 
canonical approach. Since the entropy becomes negative at low temperatures, it is necessary 
to consider a different solution for low temperatures. 

We therefore assume that < qo < . . . < qx-i < qx = 1- Then cr 2 — > for i = 
0, . . . , K — 1 and ok — ► oo in the limit p — ► oo. A simple calculation yields, 



9K-i{y) 



ln(2 cosh (3m K y) + -f3m K {Jo K ) 2 + o( e -^ J<rKmic)2/2 aT- u 



K 



(56) 



The error is exponentially small and may be safely ignored. The rest of calculation proceeds 
as before using the expansion (|5Uj) and we arrive at 



F 
N 



(3J 2 
4 

K-X 

-E 

i=0 
1 

ftrriK 



A 



1 + (P - 1) E( m i+i _ m *)<?f - Ma 1 



i=0 



+ ym 2 - ^^^(Jdx) 2 



■[rriKisech 2 (3rriK(H + Jqw)) + mj(tanh 2 (5rriK{H + Jo 7 ™))] 



(In 2 cosh /3m^ (if + J Q m)) + 0(er o , . . . , a K _ 1: a Q a 1} . . . , (T Q a K _^). (57) 



Stationarity with respect to the variational parameters gives, in the limit p — > oo, 

m = (ta.iah8m K (H + Jo m ))> (58) 

<?o = 9i = ■ ■ ■ = Qk-\ = (tanh 2 f3m K {H + J m)}, gjc = 1, (59) 
consistent with initial assumption q K = 1, and 

4 4 = (hi 2 cosh /3m^ (if + J m)) - (3m K ((H + J m) tanh 8m K (H + J m)}. (60) 

These results are the same for all if > 1, showing that no other solutions are possible beyond 
one-step replica symmetry breaking. The free energy in the limit p — > oo is given by 

— = -^~ m K + —m 2 - — (In 2 cosh (3m K (H + J m)). (61) 

The entropy is 

S (PJ) 2 1 

— — = rriK—(3 {(Jom + H) tanhPrnKUom + H))-\ (\n2 cosh (3mK(Jo m + H)) . 

Nkb 4 mjf 

(62) 

Taking into account the self-consistency equation (|5Uj) we find that the entropy vanishes 
identically. Thus this solution corresponds to the frozen phase found in the microcanonical 
approach. 

The self-consistency equations (|HHj) and (f6T)|) imply that /3rriK is independent of temper- 
ature. Since this solution is acceptable only for nix < 1, we have 

f3m K = (3 C , (63) 

where (3 C is found from the equations 
(8 I) 2 

^ J - = (In 2 cosh 8 C (H + J m)) - 8 C ((H + J m) tanh 3 C (H + J m)), (64) 

and 

m = (tanh 8 C (H + J m)). (65) 

Thus we see that f3 c corresponds precisely to the critical temperature for the transition 
to the frozen phase found in the microcanonical approach. The Parisi order parameter 



function q(x) 17( has two flat portions go = m2 an d qx — 1, with a discontinuous jump at 

x = m K = T/T c , 

q{x) = m 2 9 (I r - x ^+e(x--y (66) 
10 



The overlap distribution function 




is given by 



P(q) = T 5{q - m 2 ) + (l - £) 8(q - 1). (67) 

Thus the frozen phase is indeed a spin-glass phase with many pure states having minimal 
overlap between them and maximal self-ovelap jl^ . 



IV. PHASE DIAGRAMS 

The phase diagrams for the model defined by the Hamiltonian in the p — > oo limit 
were determined for two distributions of random fields which are often considered in the 
literature: The discrete bimodal distribution 

P{H) = 1 -5{H -H + a)+ l -5(H - H - a), (68) 



and the continuous Gaussian distribution function 



1151 



P(H) = _^ e ~(H-^?l2«\ ( 69 ) 

V27T<7 



In both cases the means and variances are Hq and a 2 , respectively. 

The Hq x T phase diagrams in the absence of ferromagnetic interactions (Jo = 0) and 
various values of the standard deviation a are shown in Fig. ^ and Fig. El for bimodal and 
Gaussian distributions, respectively. Notice that the case a = reduces to the REM in a 

are two phases: replica-symmetric paramagnetic phase (P) and a 
frozen spin-glass (SG) phase with one-step replica symmetry breaking. The transitions are 
determined by equation (}32|) for J = 0, 

(S J) 2 

- ^-^ L - Pc (H tanh (3 C H) + (In 2 cosh (3 C H) = 0. (70) 

These transitions are second order in the thermodynamic sense, but the Parisi order param- 
eter (}6l)|) changes discontinuously at the transition. For both distributions the main effect 
of the disorder a > is to increase the transition temperature for small H and depress 
it for large Hq when compared to the case a = 0. This effect is most pronounced for the 
bimodal distribution, shown in Fig. ^ For a > 0, as H$ departs from zero, the transition 
temperature initially decreases, reaches a minimum, and then starts to increase as a function 
of H . For o ^> J the minimum of the transition temperature occurs for H ~ a, as can 
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FIG. 1: The Jo = phase diagram for bimodal random field for various values of a. 

be derived from Eq. El For the gaussian case, on the contrary, the transition temperature 
increases monotonically as a function of H , as shown in Fig. |21 

The T x J phase diagrams for symmetric random-field distributions (H Q = 0) are shown 
in Figs. El to El For comparison the phase diagrams for the case a = 0, corresponding to 
the REM with ferromagnetic interactions, are shown by thin dashed lines. There are four 
phases: replica-symmetric paramagnetic (P) and ferromagnetic (F) phases, and frozen spin- 
glass (SG) and mixed (M) phases with one-step replica symmetry breaking. Unlike the SG 
phase, in the M phase there is a non-zero magnetization (m ^ 0). 

The transition from P to SG phase is determined by equation (JH2*|) for m = which is 
identical to equation (|70j). Since there is no dependence on J , it represents an horizontal 
line in the T x J phase diagram. This transition is second order in the thermodynamic 
sense, but the Parisi order parameter (jfifij) changes discontinuously at the transition. 

A second-order transition from the P to F phase can be determined by expanding the 
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3 - 



2 - 












1 

k B T/J 



FIG. 2: The Jo = phase diagram for Gaussian random field for various values of a. 
equation of state 

m = (tanli/3(# + J m)) , (71) 

in powers of the magnetization m Q|. For symmetric distribution of the random fields 
(H = 0) we find 



where 



m = am — bm — cm — ■ ■ 



a = (3J (l - (t&nh 2 pH 
b = h (3J ) 3 (l- 4 (tank 2 pH) + 3 (t&nh 4 l3H 



c = -7^ J o) 5 (2 - 17(tanh 2 /3#) + 30 (tanh 4 (3H) - 15 (tanh 6 /?i/ 



(72) 

(73) 
(74) 
(75) 



There is a second-order transition from P to F phase for a = 1 and b > 0. For 6 = there 
is a tricritical point, and for b < the transition is first order and can only be determined 
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FIG. 3: The phase diagram for symmetric bimodal random field. The results for a/J = 0.4 are 
shown by the solid curves. For comparison the results in the absence of random fields are shown 
by the dashed curves. 

numerically by equating the free energies of both phases. 

For the bimodal distribution of random fields, the conditions a = b = gives 



tanh" 1 = l - ln(2 + Vs) 



(3a = tanh" 1 — = - ln(2 + V3), (76) 



which determines the location of the tri critical point. This tricritical point occurs above the 
freezing transition ()7())1 only if the standard deviation a is greater than the threshold value 



y = ^m(2 + v/3) 



In (3 + \/3) - - (3 + \/3) In (2 + 



-1/2 



0.4584695507.... (77) 



Thus the phase diagrams for the bimodal distribution are qualitatively different depending 
on the value of standard deviation a. 

For a < cr c the phase diagrams do not differ qualitatively from the case without random 
fields, as shown in Fig. El All the transition lines are second order in the thermodynamic 
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FIG. 4: The phase diagram for symmetric bimodal random field. The results for a/ J = 0.8 are 
shown by the solid curves (second-order transition) and thick dashed curves (first-order transition) . 
For comparison the results in the absence of random fields are shown by the thin dashed curves. 

sense. However, across the F-M and SG-M transitions the Parisi order parameter (JrJrjj) 
changes discontinuously. 

For a > a c the phase diagrams changes qualitatively compared to the case without 
random fields. The second-order P-F transition line ends at a tricritical point, below which 
the transition is of first order shown by thick dashed line in Fig. This line was determined 
by equating the free energies of the two neighboring phases. The transition is of first order 
in the thermodynamic sense as well as in the discontinuity of the Parisi order parameter 
flfifij) across the transition. 

For the Gaussian distribution one always has b > when a — 1. Thus the P-F transition 
is always of second order and the phase diagrams and the nature of the transitions do not 
differ qualitatively from the case of bimodal distribution for a < a c , as shown in Fig. |SJ 
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FIG. 5: The phase diagram for symmetric Gaussian random field. The results for a/J = 1 are 
shown by the solid curves. For comparison the results in the absence of random fields are shown 
by the dashed curves. 

V. DISCUSSION 

We solved exactly the REM in a random field in both microcanonical and replica ap- 
proaches. We investigated in detail the phase diagrams for bimodal and Gaussian random- 
field distributions with mean Hq and variance a 2 . The Gaussian random fields do not change 
the phase diagrams qualitatively. The bimodal random fields, on the contrary, changes the 
Hq = phase diagram qualitatively for sufficiently large a by leading to a tricritical point 
and a first-order transition at low temperatures. The same conclusions were reached in 
the replica-symmetric study of the SK model in a Gaussian random field jsj and bimodal 
random fields 

For a ferromagnet it has been shown that a random-field distribution with a minimum 
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It is likely that 



at zero field leads to a tricritical point and first-order transition [14, 
this property remains true even when random interactions are included, although at low 
temperatures the emergence of the tricritical point may be forestalled by the spin-glass 
phase [lo| . 

Our results can shed useful light on the nature of phase diagram of more sophisticated 



spin-glass models with random fields 



9 



101 ]. One interesting feature of the phase diagram 



with tricritical point is the possibility of a first-order transition from paramagnetic to mixed 
phase. 
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